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Abstract 

The transition between photoproduction and deep-inelastic scattering is investigated 
in jet production at the HERA ep collider, using data collected by the HI experiment. 
Measurements of the differential inclusive jet cross-sections da ep /dE^ and da ep /drj*, 
where E* and rf are the transverse energy and the pseudorapidity of the jets in the 
virtual photon-proton centre of mass frame, are presented for < Q 2 < 49 GeV 2 
and 0.3 < y < 0.6. The interpretation of the results in terms of the structure of the 
virtual photon is discussed. The data are best described by QCD calculations which 
include a partonic structure of the virtual photon that evolves with Q 2 . 
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1 Introduction and Motivation 



In this paper jet production in electron-proton scattering in the transition region between pho- 
toproduction and deep inelastic scattering (DIS) is investigated. The results are interpreted in 
terms of parton densities of the virtual photon which are probed at a scale determined by the 
transverse momentum (pt) of the jets and which evolve with the virtuality of the photon (Q 2 ). 

In the photoproduction of jets Q, the photon can couple directly to a parton from the 
proton ("direct" interactions). However, the cross-section is dominated by interactions, so- 
called "resolved" processes, in which the photon fluctuates into a system of partons and one of 
these interacts with a parton out of the proton to produce the jets. This separation into direct 
and resolved processes can only be made unambiguously in leading order Q. At its simplest, 
the hadronic fluctuation of the photon may take the form of a quark-antiquark (qq) pair. More 
complicated structure is built up through QCD interactions. In addition to this point-like 
"anomalous" component ||, the photon can also acquire a more conventional hadronic structure, 
often modelled as a fluctuation into a vector meson (vector dominance model, VDM). The cross- 
section for jet production can be expressed as a convolution of universal parton densities of the 
proton and the photon together with hard parton-parton scattering cross-sections. The evolution 
of the photon parton densities with the scale at which they are probed can be calculated in 
perturbative QCD and has been extensively measured in two-photon interactions and recently 
at HERA §. 

In contrast, it is usual to consider that the only contribution to jet production in DIS is 
from direct interactions with the partons in the proton probed by a structureless photon at the 
scale Q 2 . However, in the small region of phase space where high pt jets are produced with p 2 
much larger than Q 2 , it is possible that the jet production may be most easily understood in 
terms of the partonic structure of the virtual photon together with that of the protonJf|j7], |8|, ^] . 
Parton densities within the virtual photon are expected to be suppressed [@, 1C , [0], |l2| with 



increasing Q 2 until direct processes dominate at Q 2 ~ p\. Measurements of the virtual photon 
structure in two-photon interactions require detection of both scattered leptons at non-zero 



scattering angles. Only one such measurement has previously been made [14]. The extensive 
Q 2 range together with the large centre-of-mass energy available at HERA enables a detailed 
study of the Q 2 evolution of photon structure. 

After a description of the data used in this analysis, different models are introduced which are 
intended to describe the photoproduction and deep inelastic scattering regimes. The inclusive 
differential jet cross-sections are then presented as a function of jet transverse energy and rapidity 
for photon virtualities in the range < Q 2 < 49 GeV 2 . Finally, a measurement of the energy 
flow in the direction of the photon is shown as a function of Q 2 to verify the existence of a 
photon remnant. 



2 The HI Detector 

A full description of the HI detector can be found elsewhere [j^] and only those components 
relevant for this analysis are described here. 

The coordinate system used has the nominal interaction point as the origin and the incident 
proton beam defining the +z direction. The polar angle is defined with respect to the proton 
direction, and the pseudorapidity is given by rj = — lntan(0/2). 

A finely grained Liquid Argon calorimeter [|l6[ covers the range in polar angle 4° < < 154°, 
with full azimuthal acceptance. It consists of an electromagnetic section with lead absorbers, 



20-30 radiation lengths in depth, and a hadronic section with steel absorbers. The total depth 
of the calorimeter ranges from 4.5 to 8 hadronic interaction lengths. The energy resolution is 
a{E)/E ~ Q.W/yE for electrons and a(E)/E ~ 0.5/\/E~ for pions (E in GeV), as measured in 
test beams [17]. The absolute energy scale is known to a precision of 3% for electrons and 4% 



for hadrons. 

A series of interleaved drift and multiwire proportional chambers surround the interaction 
point, enabling the reconstruction of charged particles in the range 7° < 9 < 165°, and the 
determination of the event vertex. A uniform axial magnetic field of 1.15 T is provided by a 
superconducting coil which surrounds the calorimeter. 



For 1994 data taking, the polar region 151° < 6 < 176° was covered by the BEMC |18[ , a 
lead/scintillator electromagnetic calorimeter with a depth of 21.7 radiation lengths. The res- 
olution was given by 0.10 /y/E (E in GeV) and the absolute electromagnetic energy scale was 
known to a precision of about 1%. In 1995, the BEMC was replaced by the SPACAL IS], 



a lead/scintillating fiber calorimeter with both an electromagnetic and hadronic section cover- 
ing the range 153° < 9 < 177.8°. The energy resolution for the electromagnetic section has 
been determined as 7.5%/\/^(GeV)©2.5% and the absolute energy scale uncertainty is 1% at 
27.5 GeV and 3% at 7 GeV ^0|. The hadronic energy scale uncertainty of the measurement in 
the SPACAL is presently about 10%. Both calorimeter sections have a time resolution better 
than 1 ns, enabling the reduction of proton beam induced background events. The BEMC and 
the SPACAL were used both to trigger on and to measure the scattered lepton in DIS processes 
for 0.65 < Q 2 < 49 GeV 2 . 

In 1994, the backward proportional chamber (BPC) was located in front of the BEMC. In 



1995, this was replaced by a four module drift chamber, the BDC [21], in front of the SPACAL. 
The polar angle of the scattered lepton was determined using the event vertex and information 
from these backward tracking chambers. 



The luminosity system consists of two crystal calorimeters with resolution a(E)/E = 0.1/ v E 
(E in GeV). The electron tagger is located at z = —33m and the photon detector at z = —103m. 
The electron tagger accepts electrons with an energy of between 0.2 and 0.8 of the incident beam 
energy, and with scattering angles 9' < 5mrad [6' = tt — 9). 



3 Data Samples and Event Selection 

The analysis is based on data taken by the HI experiment in 1994 and 1995. Three data samples 
were used, each restricted to a region of Q 2 with good acceptance: 

• Q 2 < 10~ 2 GeV 2 : The photoproduction sample in which the scattered electron is detected 
in the electron tagger of the luminosity system. The events were triggered by demanding 
an energy deposit in the electron tagger with E > 4 GeV in coincidence with at least one 
track pointing to the vertex region {pt > 500 MeV). More details of the trigger conditions 



can be found in [22]. The data sample used was a sub-sample of that collected in 1994 



and corresponds to an integrated luminosity of 210 nb . 

0.65 < Q 2 < 20 GeV 2 : The 1995 shifted vertex data sample, corresponding to 120 nb -1 
collected in a special run in which the mean position of the interaction point was shifted 
by 70 cm in the +z direction, enabling positron detection in the SPACAL down to angles 
of 178.5°. The events were triggered by requiring a cluster of more than 5 GeV in the 
SPACAL and timing consistent with an ep bunch crossing. The most energetic cluster in 
the electromagnetic section of the SPACAL was taken as the electron candidate 2C]. 



9 < Q 2 < 49GeV 2 : The standard 1994 data sample, corresponding to an integrated 
luminosity of 2 pb" 1 . The events were triggered by requiring a cluster of more than 
4GeV in the BEMC, and the scattered positron was taken as the most energetic BEMC 
cluster 125 



The error on the luminosity determination was 1.5(3)% for data taken in 1994(1995). For 
all three samples, the z position of the interaction vertex was required to be within 30 cm of 
the nominal position. In addition, to enable comparisons between the DIS and photoproduction 
data, for all data samples the inelasticity variable y was restricted to 0.3 < y < 0.6, the range 
where the acceptance of the electron tagger is well-understood. 

For events with 0.65 < Q 2 < 49 GeV 2 it was also required that the — °f an the 

reconstructed calorimeter clusters, which should be equal to twice the electron beam energy for 
a DIS event, was greater than 45 GeV. Monte Carlo studies showed that the background from 
real photoproduction events where hadronic activity in the backward region fakes a scattered 
lepton was reduced to less than 3% in the selected kinematic region. 

The event kinematics were reconstructed from the scattered electron 4- vector. Jets were 
reconstructed in the 7*p centre of mass frame using a hp clustering algorithm |?4j] . The merging 
procedure is based on the quantity yki which is evaluated for each pair of clusters: 



2(1 - cos6 ki) 2 2 

Vki = ^2 ™ m ( E k, E i) I 1 ) 



where Ek and E{ are the energies of the clusters and 9ki is the angle between them. In addition, to 
enable the association of particles to either a photon or proton remnant, two infinite momentum 
pseudo-particles along ±z are included in the clustering procedure, but excluded from the final 
jets. Particles are combined by the addition of their 4-vectors when y^i < 1. Thus E cu t sets the 
scale for the jet resolution and separates the hard jets from the beam remnants. In this analysis, 
E cu t was chosen to be 3 GeV. 

Jets were accepted with a transverse energy E\ > 4 GeV and a pseudorapidity in the range 
—2.5 < rj* < —0.5, where E\ and rf are calculated in the j*p frame with positive rj* corre- 
sponding to the incident proton direction. For the photoproduction sample, the E\ threshold 
was raised to 5 GeV to reduce the influence of multiple parton interactions. This selection en- 
sures that the energy flow around the jet axis is well described by the Monte Carlo models used 
for the acceptance corrections. 



4 QCD Motivated Calculations 

For acceptance corrections and comparisons with the measured jet cross-sections, several event 
generators were used. 



The PHOJET 1.03 [25] generator simulates all relevant components of the total photoproduc- 
tion cross-section. It is based on the two-component Dual Parton Model (2f| and incorporates 
both hard and multiple soft interactions. The photon flux is calculated using the Weizsacker- 



Williams approximation [£7], 28 1 and the hard processes are calculated using leading order QCD 
matrix elements. Final state QCD radiation and fragmentation effects are implemented us- 
ing the string fragmentation model JETSET 7.4 |[29|| . In this analysis, PHOJET was used to 
simulate quasi-real photon-proton interactions. 



DIS events were modelled using the LEPTO 6.5 @ and ARIADNE 4.08 |3l| programs. 
LEPTO includes the first order QCD matrix elements and uses leading-log parton showers to 



model higher order radiation. The ARIADNE generator uses the Colour Dipole model |3^] to 
simulate QCD radiation to all orders. A feature of this model is that the hard subprocess need 
not be generated at the photon vertex, and this can be regarded as generating "resolved-like" 
events. For both models, hadronisation is performed using JETSET. 



HERWIG 5.9 [33 1 was used to model direct and resolved photon processes. The emission 
of the photon from the incident electron is generated according to the equivalent photon ap- 
proximation The parton-parton interactions are simulated according to leading order QCD 
calculations, and a parton shower model which effectively includes interference effects between 
the initial and final state showers (colour coherence) is implemented |34j|. The factorisation 
scales were set according to the transverse momentum of the scattered partons, with a cut-off 
at p 1 t nm = 1.5 GeV. A cluster model is used for hadronisation. 

HERWIG includes the option of additional interactions of the beam remnant in the phe- 
nomenological soft underlying event model. A reasonable description of the jet profiles and the 
jet rates observed in the data was obtained with no soft underlying event for Q 2 > 0.65 GeV 2 
and with a soft underlying event in 15% of the resolved interactions at Q 2 = OGeV 2 . These 
values were used throughout the subsequent analysis. 

The RAPGAP Monte Carlo (3f|] originally developed to simulate diffractive processes, also 
includes modeling of deep-inelastic and all relevant resolved photon processes. DIS processes are 
simulated using leading order QCD matrix elements with a p 2 cut-off scheme for the light quarks 
and the full matrix element for heavy quarks. For resolved photon processes, the equivalent 
photon approximation is used to model the flux of virtual photons. Parton-parton interactions 
are calculated from on-shell matrix elements supplemented by initial and final state parton 
showers. 

Both HERWIG and RAPGAP include models for the evolution of the photon parton densities 
with Q 2 . Three approaches are considered. The first assumes no Q 2 dependence of the parton 
densities. The second uses the Drees-Godbole parameterization [O] of virtual photon structure, 



following an analysis of Borzumati and Schuler [13], in which the quark densities in a 

photon of virtuality Q 2 probed at a scale p 2 are related to those of a real photon / g u by: 



f qh *(x,p 2 ,Q 2 ) = f qh (x,p 2 )L(p 2 ,Q 2 ,Lo) (2) 

= ( 3 ) 



An analogous relation exists for the gluon density, with L replaced by 1? . The parameter uj is 
the value of Q 2 above which the suppression becomes significant. We use a value of uj 2 = 1 GeV 2 
and the GRV-G HO (DIS) |3(| parameterizations for the unsuppressed photon parton densities. 
Throughout this paper we refer to this as the DG model. The third approach is to use the photon 
parton densities from Schuler and Sjostrand |i7j (SaS) which are valid for Q 2 > OGeV 2 . In this 
scheme, the photon parton densities are decomposed into a direct, a VDM and a perturbative 
anomalous component. We will show comparisons with the SaS-2D parameterisation in the DIS 
scheme using the form for the Q 2 suppression recommended by the authors. 



All models were used with the GRV94 HO (DIS) []38| parton densities for the proton, which 
give a good description of the measured F2 for Q 2 > 1 GeV 2 pi] . 



5 Inclusive Jet Cross-Section 
5.1 Determination 



The distributions of the jet transverse energy (E£) and pseudorapidity (77*) in the 7*p centre of 
mass frame were corrected bin-by-bin for detector effects using generated events passed through 



a simulation of the HI detector based on the GEANT program [p9|1 . The bin sizes were chosen 
to keep the effects of finite resolution and bin-to-bin migration small. For the photoproduction 
data, correction factors were determined from the HERWIG DG model, which gives a good 
description of the data. The model dependence was estimated by comparison with the values 
obtained from PHOJET, which also gives a good description of the jets observed in the data. 
This model dependence is one of the largest contributions to the systematic error. For the 
DIS data, correction factors in the range 0.65 < Q 2 < 20GeV 2 were determined from the 
HERWIG DG model, and for Q 2 > 20 GeV 2 from the HERWIG direct model, both of which 
give a satisfactory description of the observed jets in these kinematic regimes. LEPTO was used 
to estimate the model dependence of the correction factors, which is again one of the dominant 
sources of systematic error. 

The other large source of systematic error arises from the uncertainty in the knowledge of the 
hadronic energy scale of the Liquid Argon calorimeter. This has two contributions; a possible 
3% variation in the energy scale between different calorimeter modules, which is included in 
the point-to-point error, and a 4% uncertainty in the overall energy scale, which affects the 
normalisation of the jet cross-sections. 

Further sources of systematic error include a 1(2)% uncertainty in the electromagnetic energy 
scale of the BEMC(SPACAL) and a 1 mrad uncertainty in the polar angle of the scattered elec- 
tron. For the photoproduction data, the uncertainty in the acceptance and energy calibration of 
the electron tagger was included. A 20(10)% uncertainty in the knowledge of the hadronic en- 
ergy scale of the BEMC(SPACAL) is also considered. The 1.5(3)% uncertainty in the luminosity 
determination in 1994(1995) affects the overall normalisation of the jet cross-sections. 

The effect of radiative corrections in DIS events has been studied using the HERACLES 
program [^] which includes complete first order radiative corrections and the emission of real 
bremsstrahlung photons for the electroweak interaction. The effect is 20-30% for jets with E\ of 
4 GeV, decreases with increasing E% and is negligible for Ef > 7 GeV. It does not significantly 
influence the conclusions and the data are not corrected for this effect. 

The corrected cross-sections obtained from the 1995 shifted vertex data and from the 1994 
data are in good agreement in the region 9 < Q 2 < 20 GeV 2 where the data samples overlap. 

5.2 Results 

Figure p] shows the inclusive ep jet cross-section da ev jdE\ for 0.3 < y < 0.6 and jets with 
—2.5 < rf < —0.5, and the values are listed in table The data are compared with the 
prediction from the HERWIG DG model, which includes a resolved component of the virtual 
photon. This is able to give a good description of the data except for jets in the lowest E% range 
when 9 < Q 2 < 49 GeV 2 . Also shown is the direct contribution to this model which accounts for 
an increasing fraction of the total prediction as Q 2 increases, but which alone cannot describe 
the measured jet cross-sections. 

The jet cross-section da ep /drf for jets with E\ > 5 GeV is shown in figure || and the values 
are listed in table ||. The data are compared to the HERWIG DG model and the direct photon 
contribution to this model. The direct photon processes alone significantly underestimate the 
jet cross-section at low Q 2 , but the data are described by HERWIG if the resolved photon 
component is included and suppressed with increasing Q 2 according to the DG model. The 
relative contribution to the jet cross-section from resolved photon processes increases towards 
the proton (+17*) direction. 



Q 2 (GeV 2 ) 


Flux factor 


< 10~ 2 


1.16 x 10- 2 


0.65 < Q 2 < 1.2 


6.56 x 10" 4 


1.2 < Q 2 < 2.6 


8.27 x 10~ 4 


2.6 < Q 2 < 4 


4.61 x 10~ 4 


4 < Q 2 < 9 


8.68 x 10" 4 


9 < g 2 < 20 


8.54 x 10~ 4 


20 < Q 2 < 25 


2.39 x 10~ 4 


25 < Q 2 < 36 


3.90 x 10" 4 


36 < Q 2 < 49 


3.30 x 10" 4 



Table 1: The flux factors used for the calculation of the 7*2? jet cross-sections (see equation ||). 

In order to study in more detail the Q 2 evolution of the photon parton densities, we factor 
out the Q 2 dependence contained in the flux of photons and calculate a j*p jet cross-section in 
each Q 2 range using: 

_ O'ep^e+jet+X , 
a~(*p-+jet + x — p (4J 



We use the Weizsacker- Williams approximation [27], |28| to calculate the flux of photons, F~\ e : 



Vmax f'^max 



F l{e = / dy I dQ 2 f lle (y,Q 2 ) (5) 

" Vmin 



2 



with 



The flux is integrated over 0.3 < y < 0.6 and Q^ax an d Qmin are the upper and lower edges of 
the Q 2 range. For photoproduction: 

1 - y 

where m e is the electron mass. 

This factorisation of the cross-section remains a reasonable approximation for high p t jet 
production if p 2 3> Q 2 0, a condition which is satisfied by the majority of our dataQ The 
numerical values of the flux factors used in each Q 2 range are listed in table |l]. 

The Q 2 dependence of the inclusive 7*7? jet cross-section at fixed jet E% is shown in figure ||. 
For El < lOGeV there is a significant decrease of the jet cross-section with increasing Q 2 . 
Also shown for comparison are the predictions from LEPTO and ARIADNE. We expect such 
DIS models to be valid when Q 2 > E£ 2 , where the photon cannot be resolved. This region 
corresponds to the 2-3 highest Q 2 ranges in figures ||a and ||b. It can be seen that both models 
give an adequate description of the data in these ranges. However, neither model can describe 
the data when Q 2 < E\ 2 and the virtual photon can be resolved. Also the predictions of these 
models differ significantly in this region. Although ARIADNE predicts a jet cross-section which 
decreases with increasing Q 2 in a similar manner to the data and is able to describe the data 
for Q 2 > 4GeV 2 , it is unable to describe the data at all E\ and all Q 2 . The prediction from 

1 We continue to use (0) as a formal definition of the 7*p cross-section for all pi and Q 2 . 



ARIADNE is sensitive to the parameters which limit the phase space for QCD radiation^ and 
therefore to the fraction of "resolved-like" events produced by the model, but we found no choice 
of these parameters which enabled the model to describe the data. 

The same data compared to a series of models which include a partonic structure for the 
virtual photon, as implemented in the HERWIG and RAPGAP generators, are shown in figure ||. 
In each case, the sum of the direct and resolved contributions is shown. The dot-dashed curve 
shows the prediction from HERWIG assuming the GRV-G HO structure function for the photon 
with no Q 2 suppression. The resulting j*p jet cross-section is almost independent of Q 2 , in 
contrast to the data which show, except for the highest E* jets, a significant decrease of the 
jet cross-section with Q 2 . The jet cross-section predicted by this model at Q 2 = OGeV 2 is 
slightly larger than that predicted for Q 2 > GeV 2 because a soft underlying event is included 
at Q 2 = OGeV 2 . Also shown are the predictions from HERWIG and from RAPGAP using the 
DG model for the virtual photon structure. The models are in good agreement with each other, 
and describe the data well except for jets with 4 < E% < 5 GeV when 9 < Q 2 < 49 GeV 2 , where 
they underestimate the measured cross-section. We note that in the DG model the photon 
parton density functions all vanish for Q 2 > p 2 , which is approximately the case in this region. 
The data are best described by the RAPGAP model using the SaS-2D parameterization of the 
virtual photon. In contrast to the DG model, the photon parton densities do not vanish when 
Q 2 > p 2 in this parameterization. 

The inclusive jet cross-section can therefore be understood if a partonic structure is ascribed 
to the virtual photon. Moreover, the observed Q 2 evolution of the jet cross-sections can be ex- 
plained by the suppression of the parton densities with increasing photon virtuality as predicted 
by QCD inspired models. For Q 2 > p 2 , the photon is effectively structureless. 

6 Photon Remnant 

The jet algorithm used in this analysis assigns particles to a photon remnant. The fraction of 
the incident photon's energy which is reconstructed in the photon remnant jet is given by: 

E* - E'* 

where E* is the energy of a particle assigned to the photon remnant and E* and E' e * are the 
energies of the incident and scattered electron respectively in the j*p frame. 

Figure |H| shows the uncorrected distribution of / in the data, as a function of Q 2 , for events 
with at least one jet with E* > 5 GeV and -2.5 < if < -0.5. At Q 2 = GeV 2 , where resolved 
photon processes dominate the cross-section, most of the events with jets also contain a photon 
remnant with a significant fraction of the incident photon's energy. Conversely, at the highest Q 2 
values, where the direct processes dominate, / is peaked at zero. Also shown are the predictions 
from the HERWIG DG model and from LEPTO after detector simulation. It can be seen that 
the distribution of / from LEPTO is peaked at zero for all Q 2 , as expected for a model which 
includes only direct processes. At low Q 2 , the data agree with the HERWIG DG model, and at 
high Q 2 they agree with the LEPTO prediction. 

The evolution of / with Q 2 is consistent with the picture of a resolved photon contribution 
which is suppressed with increasing virtuality. 

2 PARA(10) and PARA(15). We use PARA(10)=1.5 and PARA(15)=0.5 by default. 



7 Conclusions 



The inclusive ep jet cross-sections da ep /dE^ and da ep /dr/* have been measured in the kinematic 
range 0.3 < y < 0.6 and < Q 2 < 49 GeV 2 . 

Models in which the photon only couples directly to the partons of the hard scattering 
process fail to describe the data in the region E\ 2 > Q 2 , where the virtual photon can be 
resolved. Models which include a resolved component of the photon suppressed with Q 2 are in 
good agreement with the data. The best description of the data was obtained with a model 
which includes direct, VDM and perturbative contributions to the virtual photon structure. 

It has been established that the energy assigned to the photon remnant in events with high 
E\ jets is on average large at low Q 2 and decreases with increasing Q 2 , consistent with the 
picture of a resolved photon component which is suppressed with its increasing virtuality. 
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Figure 1: The differential jet cross-section da ep /dE% for jets with —2.5 < if < —0.5 and 
0.3 < y < 0.6. The inner error bars indicate the statistical errors, the total error bars show 
the quadratic sum of the statistical and systematic errors and the shaded band represents the 
correlated error from the uncertainty in the Liquid Argon energy scale. Not shown is the error 
from the uncertainty in the luminosity determination which leads to a 3% normalisation error 
for the data with 0.65 < Q 2 < 9 GeV 2 and a 1.5% normalisation error elsewhere. The data are 
compared to the HERWIG DG model (solid line) and to the direct contribution to this model 
(dashed line). 
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Figure 2: The differential jet cross-section da ep /dr]* for jets with E% > 5 GeV 2 and 0.3 < y < 0.6. 
The incident proton direction is to the right. The inner error bars indicate the statistical errors, 
the total error bars show the quadratic sum of the statistical and systematic errors and the 
shaded band represents the correlated error from the uncertainty in the Liquid Argon energy 
scale. Not shown is the error from the uncertainty in the luminosity determination which leads 
to a 3% normalisation error for the data with 0.65 < Q 2 < 9 GeV 2 and a 1.5% normalisation 
error elsewhere. The data are compared to the HERWIG DG model (solid line) and to the direct 
contribution to this model (dashed line). 
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Figure 3: The inclusive j*p jet cross-section <r 7 * p (Q 2 ) for jets with —2.5 < 7]* < —0.5 and 
0.3 < y < 0.6. The inner error bars indicate the statistical errors and the total error bars show 
the quadratic sum of the statistical and systematic errors. Not shown are the normalisation 
error from the uncertainty in the Liquid Argon energy scale, which is 15% at low E% and 
increases to 25% at high E%, and the normalisation error from the uncertainty in the luminosity 
determination, which is 3% for the data with 0.65 < Q 2 < 9 GeV 2 and 1.5% elsewhere. The 
data are compared to LEPTO (solid line) and ARIADNE (dashed line). 




Figure 4: The inclusive j*p jet cross-section <r 7 * p (Q 2 ) for jets with —2.5 < 77* < —0.5 and 
0.3 < y < 0.6. The inner error bars indicate the statistical errors and the total error bar shows 
the quadratic sum of the statistical and systematic errors. Not shown are the normalisation 
error from the uncertainty in the Liquid Argon energy scale, which is 15% at low E% and 
increases to 25% at high E%, and the normalisation error from the uncertainty in the luminosity 
determination, which is 3% for the data with 0.65 < Q 2 < 9GeV 2 and 1.5% elsewhere. The 
data are compared to HERWIG with no suppression of the photon structure function with Q 2 
(dot-dashed line), the HERWIG DG model (dashed line), the RAPGAP DG model (dotted line) 
and RAPGAP with the SaS-2D photon structure function (solid line). 
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Figure 5: The uncorrected distribution of the observed fraction of the incident photon's energy 
which is reconstructed in the photon remnant jet (/) for events with at least one jet with E* 
> 5 GeV and —2.5 < rf < —0.5. The data are compared to the HERWIG DG model (solid line) 
and to LEPTO (dashed line) after detector simulation. 



Q 2 




&a ep /&El 


8{stat) 


+5(sys) 


—8{sys) 


±5(norm) 


(GeV 2 ) 


(GeV) 


(nb/GeV) 










Q 2 < W-' 2 


5-7 


40.7 


0.5 


4.8 


4.0 


7.1 




7-10 


6.4 


0.2 


1.0 


1.0 


1.3 




10 - 20 


0.31 


0.02 


0.05 


0.06 


0.06 


65 < O 2 < 1 2 


4-5 


3.5 


0.3 


0.4 


0.6 


0.6 




5-7 


1.1 


0.1 


0.1 


0.2 


0.2 




7-10 


0.17 


0.04 


0.04 


0.09 


0.05 




10 - 20 


0.016 


0.008 


0.012 


0.009 


0.005 


1 2 < O 2 < 2 6 


4-5 


3.4 


0.3 


0.5 


0.6 


0.5 




5-7 


1.2 


0.1 


0.2 


0.2 


0.2 




7-10 


0.22 


0.04 


0.04 


0.06 


0.03 




10 - 20 


0.017 


0.007 


0.009 


0.007 


0.007 


2 6 < 2 < 4 


4-5 


1.5 


0.2 


0.4 


0.3 


0.2 




5-7 


0.48 


0.08 


0.16 


0.10 


0.13 




7-10 


0.049 


0.016 


0.023 


0.018 


0.008 




10 - 20 


5.1 x 10~ 3 


3.8 x 10" 3 


5.3 x 10~ 3 


6.9 x 10~ 3 


1.3 x 10~ 3 


4 < Q 2 < 9 


4-5 


1.6 


0.2 


0.5 


0.2 


0.3 




5-7 


0.69 


0.09 


0.26 


0.12 


0.22 




7-10 


0.13 


0.03 


0.04 


0.05 


0.01 




10 - 20 


0.024 


0.011 


0.019 


0.016 


0.006 


9 < 2 < 20 


4-5 


1.10 


0.03 


0.12 


0.16 


0.12 




5-7 


0.48 


0.01 


0.07 


0.11 


0.08 




7-10 


0.143 


0.007 


0.032 


0.037 


0.022 




10 - 20 


1.13 x 10~ 2 


0.09 x 10~ 2 


0.44 x 10~ 2 


0.42 x 10~ 2 


0.23 x 10 -2 


20 < O 2 < 25 


4-5 


0.26 


0.02 


0.06 


0.08 


0.04 




5 — 7 


0.107 


0.007 


0.028 


0.029 


0.016 




7-10 


0.033 


0.003 


0.017 


0.019 


0.007 




10 - 20 


5.7 x 10~ 3 


1.0 x 10 -3 


4.9 x 10~ 3 


60.8 x 10~ 3 


2.3 x 10~ 3 


25 < O 2 < 36 


4-5 


0.33 


0.02 


0.06 


0.07 


0.04 




5-7 


0.092 


0.005 


0.056 


0.019 


0.011 




7-10 


0.029 


0.002 


0.020 


0.009 


0.004 




10-20 


5.2 x lO- 3 


0.6 x 10- 3 


2.8 x 10~ 3 


2.3 x 10~ 3 


0.7 x 10~ 3 


36 < Q 2 < 49 


4-5 


0.19 


0.01 


0.08 


0.05 


0.03 




5-7 


0.051 


0.003 


0.042 


0.013 


0.004 




7-10 


0.032 


0.003 


0.014 


0.011 


0.006 




10 - 20 


5.4 x 10~ 3 


0.7 x 10" 3 


3.0 x 10~ 3 


2.3 x 10~ 3 


1.1 x 10" 3 



Table 2: The inclusive differential jet cross-section da ep /dE^ for jets with —2.5 < rj* < —0.5 in 
the 7*p centre of mass frame measured in the range 0.3 < y < 0.6 for nine different Q 2 ranges. 
The statistical, positive systematic, negative systematic and normalisation errors are given. In 
addition, the uncertainty in the luminosity determination leads to a 3% normalisation error for 
the data with 0.65 < Q 2 < 9 GeV 2 and a 1.5% normalisation error elsewhere. 
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— Z.O <. // <. 


9 1 


2.3 


0.4 


1.1 


0.7 


0.5 




9 1 << -n* << 
— Z.l <. (/ <- 


1 7 


1.7 


0.3 


0.4 


0.5 


0.4 




1 7 «r" -n* << 


1 "i 
— 1.0 


1.0 


0.2 


0.4 


0.3 


0.3 




— 1 3 < r?* < 


—0.9 


1.3 


0.3 


0.4 


0.6 


0.3 






— n 5 

u.o 


1.3 


0.3 


0.2 


0.8 


0.2 


2 6 < 2 < 4 


— Z.O <. // <. 


9 1 
— z. 1 


0.75 


0.21 


0.37 


0.25 


0.18 




9 1 << -n* << 
— Z.l <. (/ <- 


1 7 


0.51 


0.16 


0.23 


0.14 


0.05 




1 7 << -n* << 


1 ^ 
— 1.0 


0.49 


0.16 


0.17 


0.26 


0.15 




— 1 3 < 77* < 


—0.9 


0.62 


0.22 


0.25 


0.40 


0.16 




— D Q <" n* <r 
u.y \ '/ 


— n ^ 

u.o 


0.35 


0.13 


0.24 


0.16 


0.12 


4 < Q 2 < 9 


— Z.O <. // <. 


9 1 
— z. 1 


1.5 


0.3 


0.9 


0.5 


0.5 




9 1 << n* << 
— Z.l <. (/ <- 


1 7 


1.5 


0.4 


0.6 


0.6 


0.3 




— l.i V. 7j <- 


1 ^ 
— 1.0 


0.79 


0.21 


0.37 


0.28 


0.27 




— 1 3 < 77* < 


— 9 


0.64 


0.19 


0.39 


0.22 


0.14 




U.y \ '/ 


— n ^ 

u.o 


0.67 


0.21 


0.27 


0.21 


0.13 


9 < O 2 < 20 




9 1 
— z. 1 


1.25 


0.06 


0.25 


0.44 


0.20 




9 1 s' m* s' 


1 7 
— I./ 


0.78 


0.04 


0.18 


0.19 


0.14 




1 7 << n* << 


1 ^ 
— l.o 


0.62 


0.03 


0.14 


0.14 


0.10 




— 1 3 < T7* < 


— 9 

U. i7 


0.68 


0.04 


0.19 


0.23 


0.12 






— u.o 


0.48 


0.03 


0.14 


0.14 


0.09 


20 < O 2 < 25 




9 1 
— z. 1 


0.23 


0.02 


0.10 


0.11 


0.03 




9 1 ^ ti* 


1 7 
— I./ 


0.26 


0.03 


0.13 


0.15 


0.04 




1 7 << n* ^ 


1 ^ 
— l.o 


0.21 


0.02 


0.10 


0.11 


0.03 




— 1 3 < r?* < 


—0.9 


0.08 


0.01 


0.05 


0.05 


0.02 






— u.o 


0.11 


0.01 


0.07 


0.07 


0.02 


25 < O 2 < 36 




9 1 
— z. 1 


0.25 


0.02 


0.14 


0.07 


0.03 




9 1 << -n* << 


1 7 
— 1. / 


0.15 


0.01 


0.22 


0.05 


0.02 




1 7 ^ s' 


1 1 
— l.o 


0.16 


0.02 


0.08 


0.05 


0.02 




-1.3 < 77* < 


-0.9 


0.13 


0.01 


0.07 


0.05 


0.02 




-0.9 < 77* < 


-0.5 


0.13 


0.02 


0.07 


0.06 


0.02 


36 < Q 2 < 49 


-2.5 < 77* < 


-2.1 


0.37 


0.03 


0.15 


0.21 


0.03 




-2.1 < 77* < 


-1.7 


0.16 


0.02 


0.07 


0.06 


0.02 




-1.7 < 77* < 


-1.3 


0.067 


0.008 


0.087 


0.028 


0.010 




-1.3 < 77* < 


-0.9 


0.077 


0.010 


0.18 


0.032 


0.019 




-0.9 < 77* < 0.5 


0.09 


0.01 


0.06 


0.04 


0.02 



Table 3: The inclusive differential jet cross-section d<7 ej) /d77* for jets with E\ > 5 GeV in the 
7*p centre of mass frame measured in the range 0.3 < y < 0.6 for nine different Q 2 ranges. 
The statistical, positive systematic, negative systematic and normalisation errors are given. In 
addition, the uncertainty in the luminosity determination leads to a 3% normalisation error for 
the data with 0.65 < Q 2 < 9 GeV 2 and a 1.5% normalisation error elsewhere. 



